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The contribution of the box and crossed two-pion-exchange diagrams to proton-proton scattering at 90° „j is calcu- 
lated in the laboratory momentum range up to 12 GeV/c. Relativistic form factors related to the nucleon and pion size 
and representing the pion source distribution based on the quark structure of the hadronic core are included at each 
vertex of the pion-nucleon interaction. These form factors depend on the four-momenta of the exchanged pions and 
scattering nucleons. Feynman-diagram amplitudes calculated without form factors are checked against those derived 
from dispersion relations. In this comparison, one notices that a very short-range part of the crossed diagram, neglected 
in dispersion-relation calculations of the two-pion-exchange nucleon-nucleon potential, gives a sizable contribution. In 
the Feynman-diagram calculation with form factors the agreement with measured spin-separated cross sections, as 
well as amplitudes in the lower part of the energy range considered, is much better for pion-nucleon pseudo-vector 
"*7^ vis a vis pseudo-scalar coupling. While strengths of the box and crossed diagrams are comparable for laboratory 
momenta below 2 GeV/c, the crossed diagram dominates for larger momenta, largely due to the kinematics of the 
^ fcrossed diagram allowing a smaller momentum transfer in the nucleon center of mass. An important contribution 
^ arises from the principal-value part of the integrals which is non-zero when form factors are included. It seems that 
^ the importance of the exchange of color singlets may extend higher in energy than expected. 
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I. INTRODUCTION 

The nucleon-nucleon interaction at intermediate energy (up to 12 GeV/c laboratory momentum, PLab) has been 
the focus of a number of experimental 0,|2,|^i4|5|,||^0,IE0 and theoretical 0,0,0,13 studies. The momentum 
dependence of the spin transfer shows a very interesting behavior. The spin correlation observable, Cnn, for proton- 
proton scattering at 90° ,„ , has values near unity at low energies, decreases to a constant value of around 0.07 from 4 
GeV/c to 8 GeV/c and then increases to values around 0.5. 

While one-pion exchange has been shown to be a very important contributor to the NN interaction at low energies 
[l3 . ll5l[T^ and at higher energies at small momentum transfer [l3|, it alone predicts a constant value of 1/3 for Cnn- 
Not only is this value in disagreement with the data, but it is not expected that single pion exchange will still be 
important at this higher momentum transfer (see Fig. ^and comments below). 

Brodsky and Farrar 0| have shown that the simple quark-exchange mechanism also gives 1/3 for Cnn- No 
calculations of absolute cross sections with quark-exchange models exist to our knowledge, although predictions do 
exist 01 for the energy dependence of the cross section. These predictions for the slope are in quite good agreement, 
not only with proton-proton scattering, but with other scattering processes at high energies. 

While it is natural to think that this behavior could be an indicator of the nature of the elastic scattering process, 
this idea of the identification of a mechanism has languished for many years for lack of candidate theories. We present 
here the calculation of the contribution of two-pion exchange in this energy region. We find that this mechanism 
predicts the right size for the cross section, hence provides a candidate theory for the dominant mechanism in this 
intermediate energy region. 

The interest in two-pion exchan ge a s a contributor to the nucleon-nucleon interaction is very old, beginning just 
after the discovery of the pion [1^ |23, lU [22, 123 ■ The seminal work of Partovi and Lomon [2J| was one of the first 
to create a viable potential based on this idea. Since these works were aimed at obtaining the two-pion-exchange 
contribution to the NN potential, non-relativistic approximations were made. 

The Paris potential group [25l |2^ |27| has worked extensively with this contribution within a dispersion relation 
approach. Again, the interest was primarily to obtain a non-relativistic potential for the nucleon-nucleon interaction. 

Studies of the relation of chiral symmetry to the two-pion-exchange potential have been made recently |2^ which 
led to more recent studies by Rentmester et al. on the long range part of the two-pion-exchange potential. Ref. 
concludes that there is strong evidence for the existence and importance of two-pion exchange. 

Seeing this significant body of study at low energy, it is natural to ask about the two-pion interaction at higher 
energies. We will investigate its role by calculating the lowest order Feynman graphs 31] for this process in the range 
Piah from to 12 GeV/c. 

Before beginning the calculationSjSome discussion is needed as to the general approach. A common methodology 
has been to proceed by ranges [l^, llfl ■ In this point of view one says that the one- and two-pion exchange are valid 
(and give essentially all of the potential) beyond a certain inter-nucleon separation. In order to confront the data it is 
then necessary to construct the potential at shorter ranges, perhaps with phenomenology as in Ref. ||2^. We take the 
point of view that the pion exchanges are between quarks and hence (to some approximation) the interaction calculated 
remains valid even at short distances. There is no reason a priori to exclude pion exchange at short distances provided 
that the distribution of quarks within the nucleon is properly taken into account. Even when the relative coordinate 
between the centers is zero the range of the exchange between the quarks will be of the order of the size of the nucleon. 
In fact, there are a number of studies |3^ of baryon spectra based on the exchange of mesons between quarks within 
the interior of a single nucleon. 

Of course, we do not expect that one- and two-pion exchanges represent the entire interaction, other contributions 
are certainly present. We calculate the two-pion contribution with a view to examining to what extent the scattering 
properties can be explained with this mechanism alone. 

Since a common view is that nucleons and pions have an intrinsic finite size, the present approach includes the 
form factor resulting from this distribution in the Feynman integrals. There is a long history of treating form factors 
of complex systems in nuclear physics. It is often regarded as the amplitude for the probability of being able to 
transfer the three-momentum represented by its argument to the object while leaving it intact, hence is related to the 
probability of an elastic scattering or absorption. 

The inclusion of such a form factor is not simple, however. The full Feynman integrals arc in four dimensions, 
requiring a form factor with a dependence on all 4 variables. Many form factors deal only with three- momentum, 
although the problem of adding a fourth variable has been addressed by a number of groups (see section IIVI for a 
discussion). We approach the problem of the form factor from the point of view of invariants and analyticity with 
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FIG. 1: Values of 2|6p (see Ref. 0| and appendix A) at 90° „i., extracted as discussed in the text, compared with one-pion 
exchange and an estimated unitary limit. The dash-dot curve gives the energy dependence predicted by Brodsky and Farrar 
[l^ normalized at 6 GeV/c. 



guidance from the argument that the function should represent a boost into the center of mass frame of the nucleon 
in the on-shell limit. We will argue that a selection can be made on this dependence according to the view of the 
physical origin of the size of the system, i.e. if it arises from the interaction with virtual mesons (or more generally the 
sea) or with the "permanent" quark core constituents. We assume here that the pion exchanges take place between 
quarks contained within the core of the nucleon. 

We will treat the cases of both pure pseudo-scalar (PS) an d p ure pseudo-vector (PV) coupling even though some 
authors have found evidence for a mixture of the couplings [33, |3J, |33 . Some of the integrals involved in the P V 
calculation are divergent and, while the PS integrals converge, the resulting cross sections are orders of magnitude too 
large at the higher energies if no form factor is included. The introduction of form factors solves divergence problems 
which arise in the calculation but, more importantly, it takes into account the confinement core of the nucleon explicitly 
introducing an interaction range of this size. 

To represent the proton-proton amplitude, we use the Saclay |36j| decomposition into components called a, 6, c, d and 
e. These independent amplitudes have individual characteristics similar to common non-relativistic representations 
(see appendix A in Ref. J,7]) and hence provide a natural extension of our familiar concepts of the spin dependence 
of the amplitude at low energy. At 90° center of mass scattering angle, some simplification occurs since a = and 
h = — c so that only 3 independent amplitudes are needed. 

For data comparison we combine the values of Cjvat Q 111 with the measured cross sections at 90° 0, 0, Q 
to obtain the value of the absolute square of the Saclay h amplitude and the sum of the absolute squares of the d 
and e amplitudes (see Appendix A). These experimental measurements will constitute the primary data with which 
we will compare. Figure ^ shows the values of 2|&p compared with a calculation of the corresponding values from 
one-pion exchange with a dipole form factor with a range of A = 1.4 GeV/c (see Section Hvl for a discussion of this 
quantity). The one-pion-exchange mechanism is clearly important below 1 GeV/c but the predicted cross section 
decreases rapidly above that. 

There are other problems with the comparison with data of a simple perturbative calculation of the type that we 
make here. The results will not, in general, be unitary. At low energies it is not reasonable to compare with data 
without a unitary theory, which is the reason that a potential is normally calculated and followed by a solution of 
the Schrodinger equation or a relativistic generalization. We have made an estimate of what might be expected for a 
maximum cross section at 90° based on unitarity by setting each partial-wave S-matrix element to -1 and performing 
the partial-wave sum up to kc.m.R where R was taken as 1 fm. From this estimate we see that above about 1.5 GeV/c 
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the experimental cross section is sufficiently below this limit that unitarity can be expected to play a minor role. The 
introduction of a form factor also leads to questions about causality [s^lssLls^ which we discuss in section Hvl 

Section II gives the general method of calculation with the introduction of the pseudo-scalar and pseudo-vector 
operators for the box and crossed kinematics. Section III develops the method for the numerical treatment of the 
singularities of the propagators. Section IV introduces the form factors used and discuses their physical basis. Section 
V presents the dispersion relation calculation. Section VI gives the basic results of the study while section VII gives a 
summary of the work and states some conclusions. Appendix A gives the projection of the spin-dependent amplitudes 
onto the Saclay amplitudes, appendix B outlines the general method used to treat the singularities numerically and 
appendix C gives an interpretation of the variable used in the form factor in terms of a four-dimensional cross product. 



II. METHOD OF CALCULATION 



The differential cross section in the center of mass (cm.) is written in terms of the matrix element, as 




Here E is the energy of one proton in the cm. and m is the proton mass. Spin sums are implicit. Since we consider 
both the box and crossed diagrams, Al will be the sum of the two. We will study only the cases of pure pseudo-scalar 
or pure pseudo-vector coupling. 



A. Pseudo-scalar Coupling 

1. Box Diagram 

We may write the Feynman diagram for the matrix element with the box kinematics (see Fig. ^ as 



Mps5{ki +k2-k[-k'^) = 



{2n) 



dqdq'dpdp' [u{k[)j5(J) + m)j5u{ki)]-^ [uikl^h^if + m)j5u{k2)]2 



(p2 



ie)(p'2 



+ ie){q'^ — iJ? + ie){q'^ — iJ? + ie) 



T{ki,k2,k[,k'^q,q') 



xS{q-p + ki)5{q' -p + k[)S{-q - p' + k2)S{-q' - p' + fc^). (2) 

2 

Here g is the pseudo-scalar coupling constant with the normalization |^ = 13.75 and Tiki, k2,k'i,k'2, q, q') is the form 
factor derived from the intrinsic size of the interacting system. Since J- is assumed to have no spin dependence and has 
no poles on the real axes, it plays no direct role in the treatment of singularities or spin reduction. We suppress the 
arguments of T in the equations for the remainder of this section and for the next section. Using the first and third delta 
function to eliminate the integral over p and p' and the (one dimensional) relationship 5{x)5{y) — 25{x -f y)S(x — y), 
the remaining two delta functions can be written as 

2'^S{ki +k2~ k[ - k'2)S{2q - 2q' fci - fca + fca " K)- 

The first delta function factors out of the integral and cancels the one on the left-hand-side of Eq. 
We work in the center-of-mass system where 

ki^-k2=k, Ei=E2, k'i=-k'2=k', E[^E'2, (3) 

since the external lines are on shell and 



El — E2 — E 



(4) 



A Pseudo-scalar Coupling 



p= (S + go,k + q) 
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FIG. 2: Kinematics for the box diagram. Note that in our calculation, starting from Eq. ^ form factors are used at each 
vertex even though we have not indicated their presence here. 



by the conservation of energy. With these relations, the remaining delta function under the integral sign becomes 

2^S{2q - 2q' + 2k - 2k') ^ 5{q- q + k- k'). 

For the box, (see Fig. [5J|, we have 

p = k + q = k' + q', p'= k q= k' q', p' = p, (5) 

and 

90 = Po = E + qQ; p'g = E - qo. (6) 
Thus, the expression to be evaluated reduces to a 4-dimensional integral 

f dq[u{E,k')j5{^ + m)j5u{EM'HE,-\^')To{f + fnh5n{E,~k)f ^ 
(27r)4 y (p2_r7i2 + ze)(p'2_m2 + ze)(g2-^2^.jg)(g/2_^2 + j£) ■ 

We can reduce the spinor algebra of the operators on the two lines using the general relationship (for on-shell spinors) 



u{Ey>,v')'-f5u{Ey,v) = X 
where x is a Pauli spinor and we have used 

u{Ey,V 



^ VjEy' +m){Ey+m) 
2m 



Ev + m Ey> + m 



Ey + m f 1 



(8) 



(9) 



Here Ey = + m? and er is the Pauli matrix. 

The expression can be written as an operator in spin space in the form 



M 



PS 



dq ei(po,p)B2(Po.P')-^ 



With 



(27r)** J (p2 _ 77i2 ^_ ie)(p'2 — m? + it){q^ — p? -\- ie){q'^ — p? + it) 
Wr'(£^,k')75(;} + TO)75U,.(£;,k) = Ur'{E,\i!){~loPa + rn + 7 • p)ur(£',k). 



(10) 



(11) 



A Pseudo-scalar Coupling 
fci = (£,k) 



P={E- go,k' - q) 
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<? = ('7o,q) 



g' = (go , k - k' + q) 



fca = {E, -k) 



p' = {E~ go, -k - q) 



fc^ = {E,-k') 



FIG. 3: Crossed diagram kinematics. Form factors are to be included here as in Fig. |5| 

Oi can be written as an operator between Pauli spinors Xr' and Xr where the index r corresponds to the spin of the 
proton, 



^ / si Ci • k'cTi • k\ , , / cTi • k'<Ti • k 

2mei = -po{E + m) [ 1 + , ) +m{E + m)[l- ] ' 



[E + mY 



+ {E + m)[ 1 -f^ 



^i-pW 1 



2mei = i;+(TO -po) ^+ ("1 +Po) + cTi • pcTi • k + (Ti • kVi • p, 



2mQi = E+{m - po) - ™ (k • k' + icri • k' x k) + p • (k + k') + icri • p x (k - k'), 

where E^ — E + m. These expressions allow us to identify G and H in the operator, 

2mei = G + i(Ti • H, 



(12) 
(13) 



G = £;+(m ~ po) - ^^^k . k' + p . (k + k'), 



H 



For the lower line where 



I^k'xk-px(k-k') 



2me2 - G' + i(T2 ■ H', 



G' = £;+(m-p^)-^^^k.k'-p'.(k + k'); H' 



^^^^^k' X k + p' x (k - k') 



(14) 
(15) 
(16) 
(17) 



B Pseudo-vector Coupling 
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2. Crossed Diagram 

The expressions for the matrix element in Eq. H10|l . and for Qi and 02, for the two-pion crossed diagram are the 
same as in the case of the box except that now the relations between momenta arc those corresponding to the diagram 
in Fig. □ 

p + q' = k, p' + q= k, p + q = k', p' + q' = k', (18) 



pq=Po; 90 = 9o; qo = E-po. (19) 

Since only neutral pions can be exchanged in the case of the box diagram and both neutral and charged pions can be 
exchanged in the crossed diagram, a factor of 5 multiplies the crossed diagram result. This isospin factor is explicitly 
given in the dispersion relation approach in Eqs. (|169|l and (|170|l . 

B. Pseudo-vector Coupling 

In this case the interaction is given by j^'^'ltilsT ' d^'t'T^'^ s-nd for one nuclcon propagator, Ppv, the expectation value 
leads to 



/■2 

Ppy(p2 - ^ ie) = -l-u{k')i'-iz{i> + m)i-ir^u{k). (20) 



Using k + q= p = k' + q' and the Dirac equation {jk — m)u{k) = we can write 



f2 

Ppvip'^ - m'^ +ie) = ^u{k')jr,{^ + m){^ + m){-^ + m)-f5u{k). (21) 



By regrouping the terms 



f2 

PpviP^ - +ie) = ^u{k')j5{^ + m)[{:^ - m) + 2m][{^ - m) + 2m]j5u{k), (22) 



we find 



f2 f24m^ 
Ppvip^ - + ie) = ^u{k')-f5{p^ - m'^){^ + 3m)-f5u{k) + - — —u{k')-f5{^ + m)75u(fc). (23) 



The last term can be seen to be the pseudo-scalar expression considered in the last section for one nuclcon propagator. 
Hence, the operator can be separated into a term which corresponds to a contact term and one which is identical to 
the pseudo-scalar expression developed before. 

Ppv = ^-,u{k'),,ii> + SmHuik) + g^^J^J^l^Ipm ^ p, + p,,, (24) 



where we have used the standard correspondence = . We can write 

4^2 '^rriUr'{k')ur{k) + xl''S>i{po,p)Xr ■ (25) 

It is useful to define an operator analogous to the 8i(po, p) as before for the contact term for the upper line for the 
case of the box diagram 

Xr'Cl(PO,P)Xr = 2mUr'{k')Ur{k) 4- Xr' ©1 (PO , P)Xr • (26) 



Pc= ' 



A Contact-Contact Term (Bubble Diagram) 
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We can also define operators corresponding to the contact terms which are very similar to those defined before in Eqs. 
(UH) and lUSl). 



Gc--G + 2m{E+ - : G'r = G' + 2m ( E+ ^' ' ^ 



E+ 



(27) 



with 



E+ 



E^ 



2mCi(po,p) =Gc+ i(Ti ■ He; 2mC2(po,p) = G'c + i(T2 ■ 



(28) 



(29) 



The same generalization to the crossed diagram can be made as in the previous section and we do not repeat it here. 
We can write for either the box or crossed diagram 

9^ 1 f , C,{p^,v)C2{p'^M)T , g4 1 . Ci(po,p)e2(p;„p')-^ 



16m4(27r)4 7 ^ [ql - w^){ql - lo'^) 4m^ {2tt)'^ J ^ (p'^ - m^){ql - uj'^){ql - lj'^) 



dq- 



ei(po,p)c^2(p^,p')-^ 



M 



PS- 



(30) 



4m2 (27r)4 7 ^ (p2 _ m2)(g2 _ ^2)(q2 _ ^,2) 

The first term will be referred to as the bubble diagram and the second and third (numerically equal) as the triangle 



III. TREATMENT OF SINGULARITIES 

In principle, each of the propagators contribute two ^-functions and two principal- value integrals. However, for one 
fcrmion on shell neither pion can be on shell and for one pion on shell neither fermion can be on shell. This fact 
considerably reduces the number of terms which contribute. For numerical calculation, the poles and their order must 
be identified and separated into 5- function and principal value parts. The general technique for doing this is outlined 
in Appendix B. We cannot use the classical methods of reduction ^] because of the existence of the form factor at 
each vertex. In the following we treat the case for a general scattering angle x assume that the incident beam 
direction, parallel to k, is along the z-axis and the final proton momentum, k', lies in the z-x plane. 

A. Contact-Contact Term (Bubble Diagram) 

For the first term in Eq. (|30|l for either the crossed or box configuration we have 



Mcc = 



9^ 1 r , Ci(po,p)C2(p^,p).F 



16m4 (27r 



{ql-u'' + it){ql-^''^^it) 



(31) 



d£L \ dq dqo- 



0{qo,q)T 



16m4 (2^)4 J -"J, (g2 _ ^2 + ,,)(^2 _ ^,2 + ' 

where g = |q| and the volume factor q^ has been included in 0{qo, q). Taking large finite limits for the integrals 



(32) 



dq 



1o 



dqoO{qo,q) 



1 



(7q — cj2 + ie q^ — lu'^ + ie 



T. 



(33) 



B Contact-PS and PS-Contact Terms (Triangle Diagram) 



9 



If we consider the contribution of one pole in qg [qq — uj for example) then a counting of powers of q in the numerator 
and denominator shows that, for large g, the integrand is proportional to q so that the integral diverges as g^^j. if the 
form factor is set to a constant. The integral can now be evaluated by the method given in Appendix B. 
For the first term of the integrand of the solid angle we have 



dq 



1o 



0{qo,c\)T 

"■10— 5— — 

— w"^ + le 



(34) 



dq 



where 



dqo 



Oiqo,q)-T{qo,q) , r(go,q) 



qQ — oj^ + ie qQ - u}^ + ie 



T = 



O(go,q)-T(go,q) ^ 

uqo 2 9 

9o 



dq 



F{q,n), (35) 



and 



^^^ [0(c^,q)+0(-u;,q)] 



2uj 



[0{aj,(i)+Oi~iO,q)]T 



(36) 



T(go,q) 



■0{uj,q) 0{-uj,q) 



2uj 



2uj 



(37) 



We have included the +ie in the denominator of the external factor 1/ {lo"^ — uj'"^) but, in fact, there is no singularity in 
the combined integral so this choice is arbitrary. The remaining indicated integral is computed numerically. We now 
need to treat the singularity in the g integral in a similar manner. Since 



= q2 - q'2 = -2A:2(l -cosx) - 2q- (k -k') = -2kX{q-w) 



(38) 



with 



a;^ sin X cos ^ cos 6* — sin cos 0; and w = — fc(l — cos x)/X. 



X = x(l — cos x) — vl 
Here 9 and are the angles of q in polar coordinates and x — cos 6. We can now write 



(39) 



dq 



10^ — Lo'^ + ie 



Fiq,n) 



1 dqF{^,n) 
2kX Jq q — w + ie 



(40) 



2kX 



5— dq[F{q,n)-F{w,n)] 



q — W 



1 ITT 

Mq^nax/W - 1) + -—F{w, n). 



2kX 



2kX 



The imaginary parts will cancel when the two terms of Eq. (|33|l are combined. 



(41) 



B. Contact-PS and PS-Contact Terms (Triangle Diagram) 

In this case there is one nucleon propagator in addition to the two pion propagator and the second and third terms 
in Eq. H3U|) can be written in the form 
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4m2 {2TrY 



dn 



dq 



la 



dqQO{qo,q) 

uj'^ — oj'^ + ie J-q^"^ + ie 



1 



1 



9q — cj^ + ie qQ - Lu'- 



(42) 



C PS-PS Term 
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A count of the powers of g at a pole of qo shows that the integrand varies as 1 /q for large q so the integral diverges 
(in the absence of a form factor) logarithmically in qmax- Again there are no double poles and the first method of 
Appendix B can be used to calculate the 4 poles in q^. 

For the crossed-pion configuration we define, for the first term in the above equation 



while for the second we have, 



With the definition 



Pi^E-Ep; P2^E + Ep; P3 ^ uj- P4 = -w, 
Pi=E-Ep; P2 = E + Ep; P3 = P4 = -cj' ■ 



the integral can be written as 
9' 1 



dq 



4m2 (27r)4 J 
The last integral in the above becomes 



I0 



dqo 



G((7o,q) 



[O(go,q) - G(gQ,q)] _^ 

Utiido-P. + ^e) UUiqo-P. + ^e) 



«o G(go,q).^ 
dqo — 4 



9o 



0(P.,q)^ 



It (go -P+ ^e) U.^P^ P,) 



(43) 
(44) 

(45) 
(46) 
(47) 



OiP^,^)J' 

'-tn,^^ip^-pJ) 



In 



Qo ~ 
Qo + 



o{P^,^)^ 



(48) 



C. PS-PS Term 

The pseudo-scalar term is given by Eq. (|10(l . There are now two more factors of q in the denominator so the integral 
is convergent. For this case, i.e. when there are two nucleon poles, the box and crossed diagrams must be treated 
separately. 

1. Crossed Diagram 

For the crossed-pion configuration, where po — p'o — E ~ Qo, there is formally a second order pole when Ep = Ep>, 
so we write 

1 



ipl -El + ie){pl - El, + ie){ql - + it){ql - u'^ + it) 



1 



^ - El 



pt^El+ze p'o-Ep 



qg - UJ-' + ie qg - uj'^ + ie 



The integration can now be done with the methods of Appendix B with four poles in the gg integration. 
The first factor can be evaluated using 



Ep - Ep, 



2 12 
P P 



(k'-q)2-(q + k) 



1k|/(a;) 



(49) 



(50) 



C PS-PS Term 
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where 



f{x) — x{l + cos x) + Vl — cos sin x ^ — > {x + y/l — x'^ cos (f>) . (51) 
From this expression we see that the second pole has been transformed to the x integration. The singularity falls at 



Xo = 



sm X cos ( 



X=90° 



cos ( 



( 1 + COS xY + sin^ X cos^ 



and gives a delta- function contribution of the type S[f{x)] with 



/'(xo) = 1 + COSX 



sin^ X cos^ x=90° 



\/l + COS^ (/) 



1 + cos^ 



1 + cos X 

Thus, for the first of the four terms in Eq. I|49|l we can write 

1-1 fi^) Jo - J-,^^^ {pI - El + ^e){ql - io^ + ^e) 



9 



2k{2ny 



d(j) 



(52) 



(53) 



(54) 



9 



2fc(27r)^ 



d^ I ' dx^^ 



(55) 



where 



2fc(27r)4 



dx 



J(a;» - Q{x,(l)) _^ (5(a;» 



f{x)+te f{x)+: 



Q{x,(f>) 



/(^) 



— -J(a;o». 



f'{xo)ix - X0 + if) 

The first integral in Eq. H56|l now has no singularity and can be performed numerically. The second becomes 



dx 



J{xo, (j)) ^^f 1- Xo 



1 f{x)+ie f'ixo) J_iX~xo + ie /'(xo) Vl + ^^o 



(56) 



(57) 



(58) 



2. Box Diagram 
For the box diagram we have for the pole structure 

1 



(p2 ~El+ ie){p't, + ie){ql - u"^ + ie){ql - uj'^ + ie) 



The pion pole part can be expanded as before but the nucleon poles have a very different development since the fourth 
integration variable is not the same in the two factors but Ep = Ep> leading to a second order pole. We can write 
(with po = E + qo and p'g = E — qo) 

1 1 



[pI - {Ep - tef] \p'l - {Ep - ief] [{E + q^Y - {Ep - te)^] [{E - qof - {Ep - ze)2] 



C PS-PS Term 
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The first factor can be handled by standard methods since there is no second order pole. The second factor, however, 
does contain a second order pole when Ep = E. For this factor we write 



E + E,, 



ql-{E-Ep + lef 2{E^ - + «e) 



1 



qo - E + Ep - ie qq + E - Ep 



The integral can now be done with the methods of Appendix B with five poles in the qg integral in the following 
manner. 
We express 



-q^ — 2kqx = —q{q + 2kx). 



Including the pion poles with this factor we have 



E + Er, 



2q{q + 2kx + ie) 2kX{q - w) ' 



_qo — E + Ep ~ ie qo + E — Ep + ie 



where w was defined in Eq. (|39|l . 
We expand the factor 



and write 



gg — + ie (7o — uj'^ + ie 



1 1 



(60) 



{q + 2kx) [q — w) w + 2kx 



X 



q + 2kx w — q 



X 



w + 2kx k{l-2xX) kf{x)' 



where 



/(x) = 1 — cosx ~ 2xc(a;) and c{x) — x — sin xi/ 1 — x^ cos (f> — x cos x- 
The function f{x) has two zeros given by 



'1 ± 



x+ = ±1 



VT+^ smxcos( 
with z = 



2 1 — cosx 

The "+" under the radical goes with the "+" on the exterior and vice versa. The derivatives at the zeros are 



/'(^±) 



2ViT^ 

x± 



(61) 



(62) 



(63) 



(64) 



(65) 



Denoting the inner integrals over q (and go) Ji and J2 corresponding to the two terms in Eq. H61() we write the 
integral over the variable x as 



9 



4fc2(27r) 



dx 



Mx,(j>) - Q^{x,<l)) Qi{x, 



where 



i(x,^) = J{x) 



f{x)+ie f{x)+ie 
Ji{x+,(l)) Ji{x-,4>) 



f'{x+){x - x+) f'{x-){x - X-) 



(66) 



(67) 



C PS-PS Term 
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The second term in Eq. H()6|l can be evaluated by 



dx 
1 f{x 



■In 



x+ 



1+X+ 



with 



f'ix-) \l + x 



Ji{x,cf)) 



^^q + 2kxj_^ ^%l-{E + Ep-ieY 



(68) 



qo — E + Ep ~ ie qo + E - Ep + ie 



and 



q^ — u)^ + ie (7q — cj'-^ + ie 



2 _ ,,,/2 



0{qo,q)T 



w-q J-oo Qo - {E + Ep + iey 



(69) 



1 



qo~ E + Ep~ ie qo + E - Ep + ie 



ql- Lo"^ + ie ql - w'2 + ie 



(70) 



IV. FORM FACTORS 

The two protons are treated as finite-sized particles with their intrinsic internal structure providing an extended 
source for the exchanged pions. In this case the range of the form factor is not a cut-off parameter to be taken to 
infinity at the end of the calculation, as is often done in a local field theory, but a physically meaningful quantity. 
We consider that the true interaction of the pion is with the partons within the nucleons but we use an effective field 
theory (assumed to obey Feynman rules) to describe the composite system. 



Basic Considerations 



In general, one may expect a finite distribution of the interacting constituents of the nucleon to have two components: 
one due to the confinement range of the valence quarks (and other non-color-singlet partons) and one due to the meson 
cloud. Because of G parity, the exchange pions cannot couple directly to the pion part of the cloud (likely to be the 
most important part) and the cloud can be expected to have a larger extent (and hence be "softer" in momentum 
space). For these reasons we consider only that part of the density due to the distribution of valence quarks. This 
argument is by no means new (see e.g. Maekawa and Robilotta ^3 fo^' the case of one-pion exchange as well as Refs. 
|0)E3)E3)- Form factors of the mesonic dressing type have been extensively studied (see e.g. Ref. |4^). 

The non-relativistic form factor can be obtained directly by calculating the Fourier transform of the density (see 
e. g. Ref. [13). While more realistic forms are possible, we assume here an exponential density in the center of mass 
of the nucleon corresponding to a dipole form factor of the type 



(A^-/i^)^ 
(q2 -I- A2)2 



A2 



A2 



A2 



A2 



(71) 



if we chose to normalize the form factor to unity at the nucleon pole and we have neglected the small recoil correction, 
/i'*/(4m2), to — /i^. The spatial distribution corresponding to this form is proportional to 



-Ar 



It is well known that 



A Basic Considerations 
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Pion— Quark Distributions 

Dashed: A = l.2 GeV/c (K =0.570 fm) 
Solid: A = 1.4 GeV/c (K =0.488 fm) 
Chain: A = 1.6 GeV/c (R =0.428 fm) 
Open and full circles: Lissia et al. 
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FIG. 4: Comparison with distributions arising from the dipole form factor with those obtained from the lattice calculations of 
Lissia et al. psj by the procedures given in the text. 



there is a problem with this density since it has a finite derivative with respect to r at the origin but this presumably 
only causes significant corrections at very high momenta. The value of A = 1.4 GeV most often used in our calculations 
corresponds to an rms core radius of the exponential density of 0.49 fm. 

An exponential shape is indeed suggested by lattice calculations. Lissia et al. 45] gave spatial distributions for the 
distance between quarks, Pqqirqq), in the nucleon, and the pion and p mesons. To find the distribution about a fixed 
center, Pq{r) needed here, a conversion must be made. If we assume no correlations between quarks and make no 
center-of-mass correction (in the same manner as Lissia et al. did to compare with the MIT bag model) the conversion 
between these two densities is given by the rule that the Fourier transform of pqq is the square of that for pq. The 
distributions given in Ref. for pqq are well fit by an exponential, e""''"''', except for the smallest values of Vqq, 
leading to a dipole form factor. The square root of this form factor gives a monopole form factor with the same value 
of a. Since the pion-nucleon interaction will be governed by a convolution of the pion and nucleon distributions, or a 
product of the form factors, the result of this type of conversion is the product of two monopole form factors, or, if the 
ranges are the same, a single dipole form factor. The ranges we extract from the figures in Ref. are a — 7.0 fm~^ 
for the proton and f3 — 6.4 fm~^ for the pion. If we were to take both ranges to be 7.0 fm""'^ the result would be a 
dipole form factor with a range very nearly 1.4 GeV/c which is what we have used in most of the calculations shown. 
For two different values of the ranges the density is given by 



P^'^i^) = l^^- '- (72) 

— r 

This density is shown by the open circles in Fig. ^compared with densities resulting from dipole form factors, i.e. 
pure exponentials. The the root mean square radius, Rms, corresponding to this distribution is 0.52 fm. 

Another possibility for the conversion of the pion distribution is to take the distance to the center of mass to be 
half the distance between quarks. This choice leads to a product form factor 

' " > (73) 



Q,2 _|_ q2 \^ ^2 _|_ q2 



C Singularities in qo 
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where the /? is twice the value found in the fit to pqq. This form factor gives the density 



r(/32 



-0r (/3^ - a^)r 

P — f 

2/3 



(74) 



which is shown as the sohd dots in Fig. 0] The Rms corresponding to this density is 0.44 fm. 

Alexandrou et al. |4^ also made calculations of and commented that their distributions were very well fit with 
an exponential form. In a variety of calculations they found Rrms values of the pion-proton system in the range from 
0.55 to 0.60 fm. 

Another estimate of this radius can be obtained directly from form factors for pion-nucleon coupling. For a monopole 
form factor, a value of A of 800 MeV to 1 GeV was found |43|. Choosing an equivalent rms radius leads to values in 
the range 1.13 to 1.41 GeV/c for a dipole form. 

While the original MIT bag model used a radius near 1 fm. (for all of the constituents of the nucleon), corrections 
due to the inclusion of a pion cloud and the recoil degree of freedom made by Bolsterli and Parmentola 48] led to 
values of the rms radius near 0.45 fm. 



B. Relativistic Generalization 



One often attempts to make a generalization of Eq. (|71|l to a relativistic dependence on momentum of the type 



(A2-^2)2 



ill 



A2)2 



(75) 



I.e. 



= q^ — ■ While this substitution makes the form factor invariant, it is not the only way to achieve 
this objective. This procedure introduces additional singularities on the real axis in go which are often regarded as 
spurious. 

We suppose that the distribution of quarks is known in the rest frame of the nucleon and assume it to be spherically 
symmetric with its Fourier transform being a function of the square of the conjugate momentum, t^. This distribution 
should be boosted into the center of mass system of the two nucleons but a boost is only defined for on-shell particles. 

On general principle, to preserve Lorentz invariance, the form factor should be a function of Lorentz scalars only. 
In the case of a pion-nucleon vertex (neglecting spin) there are three four-vectors to work with, the initial and final 
nucleon momenta and the pion momentum. Due to 4-momentum conservation only two are independent. Let us 
choose the initial nucleon momenta (k) and the pion momentum (q). From them wc can construct three scalars 
q^ , /c^, and k ■ q. 

We expect the desired invariant to involve the nucleon momentum if it is to represent a boost into the nucleon rest 
frame in the on-shell limit. Working with these invariants and taking into account that there should be scaling with 
the dimensions of q, to evaluate to t^, the expression should be bilinear in q. Since k ■ q is linear in q we must use its 
square, which implies that, in order to be coherent in k as well, the only two invariants to consider are (k ■ g)^ and 
fc2g2 xhc linear combination 



(k-qf -k^q^ = m^Z^{k,q) 



(76) 



indeed reduces to m t — m q when the nucleon four- momentum, k, is taken to be (m, 0). We shall adopt the 
variable Z'^{k,q), which provides an off-shell generalization of the boost to the nucleon rest frame, to replace in 
Eq. (|71|l . For a discussion of the Lorentz tensor which represents the 4-dimensional cross product of k and q and can 
be contracted to form see appendix C. 



C. Singularities in go 

The problem of the possible unwanted poles in qo on the real axis is at least partially resolved since the condition 
for a singularity (fc • q)^ — k^q^ + nri^h? — gives as roots for go 



90 



fcok • q ± ^-{fc2[k2q2 - (k • q)2] + \<i^m'^K^ 

k2 • 



(77) 
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As long as the nucleon remains time-like (fc^ > 0) the quantity under the radical is negative so that there is no zero 
on the real axis. Since, for the exchange of two pions, at least one of the nucleons at each vertex is on shell, and there 
is no zero on the real axis for qq for the nucleon momentum on shell, the problem of unwanted additional singularities 
is resolved. 

Poles do exist in the complex plane and lead a finite value of the principal- value integral on qo (an integral over only 
poles on the real axis will lead to a principal value contribution of zero). Since, if the pole is very far from the real 
axis its contribution will become negligible, in the limit of A — > oo the principal-value part of the Qq integral vanishes 
as expected. We can see from Eq. (I77f) that in the limit of k ^ the pole also moves far from the real axis so at low 
energies the principal-value contribution tends to zero. 

Such poles in the upper half of the complex plane are often considered to lead to a violation of causality. In the 
present case we are considering a pair of finite-size particles which are described in terms of the motion of their centers. 
Hence, they can start to interact when they are a distance apart equal to twice their radii. The usual considerations 
of causality Is^ Is^ deal with two point particles and the expectation that any interaction between them can not 
propagate faster than the speed of light. This criterion does not apply to the present case since the points which 
need to be connected by the speed of light are not the centers but any co-moving points included in their structure. 
The scattered wave can start to appear before the centers of the distributions can be connected by a light ray. The 
appearance of these poles is natural and expected. Of course, their contribution tends to zero as the size of the system 
goes to zero (A ^ oo). 

Since the size of the system becomes unmeasurable at very low energy we should expect that the nucleons could 
be treated as point particles in this regime so that we expect that standard causality will be valid in the low-energy 
limit. We see that, since the pole off of the real axis gives a negligible contribution in this region being very far from 
the real axis this expectation is realized. 



D. Properties of Z 

The function Z^ has the property 

Z^{p±q,q) = Z^{p,q) or, in particular Z^(k',q) = Z^{k,q). (78) 

Since Z^ is independent of the use of the initial or final nucleon momentum, the vertex function is a property of the 
vertex itself and not of the individual four-momenta. If the three 4-momenta which are connected to the vertex are 
Pi, p2 and P3 with a conservation law between them, say p^ — ±pi ± p2, then any two of the momenta can be used 
for the evaluation of the vertex function, i.e. 

Z^ip,,P2) = Z2(pi,P3) = Z^{p3,P2) = ... (79) 

An equivalent way to express this function is 

^2 ^ [pI? + [pI? + [pI? - '^{pIpI+pIpI+pIpD 

In a more general view, there are two invariants available; [k ■ q)"^ jm?, which evaluates to q^ in the rest frame of 
the nucleon and '''\^2 — — which evaluates to in this rest frame. One could choose any combination of q'^ and 
g§ desired for the dependence in the rest frame. However, only [{k ■ q)^ — k'^q^]/w? (and q^ by dint of containing no 
reference to the nucleon momentum at all) are independent of which nucleon momentum (fc or k') is used. 

The use of the form factor obtained from the Fourier transform of the density in the nucleon center of mass without 
any dependence on go corresponds to an interaction which is instantaneous in time. This is perhaps to be expected 
since the valence quarks are always present for interaction, whereas the interaction through the intermediate step of 
the emission of a meson requires a time ?i/Mass for the meson to be formed and propagate. 

Other authors [i^ Isof have used relativistic generalizations of the form factor. In particular Ramalho, Arriaga 
and Pena 49] have argued that the first type of form factor dependence, Eq. (|75|l corresponds to the dressing of the 
nucleon by mesons of mass A and suggest that the variable q^ should be replaced by 

^ - q\ (81) 



F Application to the present case. 
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where P is the center-of-mass momentum of the two nucleons. The value of evaluates to in the center of mass 
system so that there is no dependence on qo if one calculates in that reference frame, hence no possible singularity in 
that variable. This may be a useful form but the underlying basis in physics is not clear since the form factor should 
be a property of the pion-nucleon vertex and not the entire system. How a pion interacts with a single nucleon should 
not depend on the momentum of a second nucleon. 



E. Conditions on the Form Factor 



Reviewing the previous paragraphs, we may summarize the desired properties of an off-shell form factor arising from 
an extended source. 

1. There should be no poles of qo on the real axis. 

2. For the case of the nucleon on shell and at rest, the form factor should be independent of qo so that the interaction 
is instantaneous. 

3. As the nucleon energy approaches zero, poles in qo should tend to ±ioo. This is necessary since, in this limit, the 
interaction could be regarded as being among point particles and, by causality, there can be no contribution from 
poles off of the real axis. While the arguments of the previous point and this one arise from different physical 

bases, they are mathematically related. 

For a typical example consider a form factor of the type 

\2 

(82) 



aqo + bqo + + ' 
where a and b are real functions of q, k and ko- The poles in qo will be given by 



-6± V62-4a(q2 + A2) 
qo = ^ . (83) 

In order that g'o ^ ±ioo as |k| — > both a and b must tend to zero in this limit. This, in turn, means that the 
form factor becomes independent of go in the limit |k| ^0. 

4. In the nucleon rest frame the form factor (which is a function of |q|^ only) is the Fourier transform of the source 
density, a positive definite function. 



F. Application to the present case. 



For the application of the form factor in the present calculation we can always choose the nucleon to be one of the 
external lines and hence on shell with energy E. If we define 



/(£;,k,go,q) = f{k,q) = 



A2-^2 



H 2 



it3)l -q2 + A2 

a product of four of these factors will appear, one for each vertex. For the box diagram we have 



(84) 



^ = f{E, k, qo, q)f{E, -k, go, cDf{E, k', go, ci')f{E, -k', go, q'), (85) 

while for the crossed diagram the product 

= f{E, k, go, ci')f{E, -k, go, cOf{E, k', go, q)/(iS, -k', go, q') (86) 

enters. Since these functions are scalars and do not have any poles on the real axis, the manipulation of the spin 
observables and the treatment of the singularities is not affected by their presence. 
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V. DISPERSION RELATION APPROACH 

In order to cross-check our Feynman diagram results we have calculated the two-pion exchange in the dispersion 
relation approach following Ref. ^5l|. The matrix element M of Eq. (1) in section II can be written as, 

5 

M^Yl [^Pt ^ + ^ ^1 ■ "^2] P^■ (87) 

Here, w, t, i are the Mandelstam invariants satisfying w + t-\-i = Arn^, '''1(2) the isospin Pauli matrices for the nucleon 
1(2) and Pi the perturbative invariants. One has, 

W={ki+k2f, t={ki-k[f, t={ki-k'2f (88) 

and 

Pi - lll^ 

P2 = -[(7^-A^)l2+(7^-P)li], 

P3 = (7' • N) (72 . P) , (89) 

P4 = 7^•7^ 

with 

N=l{k2 + k',), P=l(k,+k[). (90) 

Assuming the Mandelstam representation the two-pion exchange invariant amplitudes pf {w,t,t) of Eq. H87|l 
satisfy the following double dispersion relations, 

±. ^ 1 /■+°° dt' r+°° dw' ±, , 



" " J 4/^2 t' — t — le 74^2 w' — w — le 

where the yf{w, t) are the double spectral functions. In Eq. H91I) the first double integrals represent the contributions of 
the two-pion-exchange box diagram (see Fig. [21) . The second double integrals give the contribution of the crossed-pion 
diagram (Fig. They are obtained from the contributions of the box diagram by the substitution ki <-!■ —k'^ which 
is equivalent to that oiw^i. This substitution changes the sign of the invariants Pi , P3 , P5 while P2 and P4 are not 
modified. It also changes the isospin dependence. These changes result in the signs given here. We shall show below 
how one can calculate the two-pion-exchange double spectral functions yf{w,t). Use of the following decomposition 



t — t' — le w' ^ [Ann? ~ w ~ t) — le t — t' — le w' — {Am? — w — t') — le 

1 1 
-I . (92) 

w' — (4m2 ^ w ^ t) — It t' — {Am? — w — w') — le 

allows us to write Eq. as 
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1 df' 

pf {w, t,t) = - [pf {w, t') T (-l)Vf (4m2 ^w~t', t')] 

p±Ki)=Im,p±Ki) = - / yf{w',t), (94) 

TT J w' — w — le 

1 Z"^"" di' 

af{w'Am'-w~w')^- ^ yt{w',t'). (95) 

^ ' TT 74^2 r — (4to^ — w — w') — le 

Since the term with af (w' , Am? — w — w') in Eq. (|93|l gives rise to a rather short range force with mass exchanges 
larger than 2m it was neglected in the Paris Potential [25L l27l | . In Eq. (|94|l the notation Imf comes from the use 
of the relation 



with 



and 



1 P 



= 1- nT6(t' - t), 

t'-t~ie t'-t ^ 

where P indicates the principal value part of the integral. 

= iVi/2, k') J9 = (go - Vi/2, q - k') k[ = {Vt/2, -k') 

A A 

Q = {qo, q) I \q' = {Vt- go, -q) 

A A 

k2 = iVi/2, k) p' = {Vt/2 - go, k - q) k'^ = (v^/2, -k) 



(96) 



FIG. 5: Box diagram kinematics in the t-channel cm.. 
The unitarity condition in the t-channel for the NN 2tt ^ NN process leads to 

EJpf (^«,t) T (-l)>f (4m2 -yo- t,t)]P, = 



1287r2 V t 



1 t-Ap. 



J dQQ [Af (si,i,si)li+7i •QB*^ (si,t,si)] [A± {s2,t,S2)l^ +7^ ■ClB^is2,t,S2) 



(97) 



where the four- momenta of the on-shell exchanged pions, in their t-channel center of mass (see the kinematics given 
in Fig. inj, are, 



2 ' ' " V 2 
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The amplitudes A(_B)^2) (si(2)i^i ■^1(2)) the usual invariants for the reaction NN tttt [s^- The Mandelstam 
invariants Si(2): t and Si(2), satisfying, Si(2) +t + Si(2) — 2m? + 2y^^, are 

si = [ki-q'f = m2+^2_^_2P-Q, 
t = {k,-k[f - 4(Q2+^2)^ (99) 

si = (fci + g)^ = m2 + ^2 _ ^ _^ 2P • Q 

and 

S2 - (fe - = to" + - ^ - 2N • Q, S2 - (fc2 + - to" + ^2 _ ^ ^ 2N • Q, (100) 
with, q2 — g'2 _ ^2 _ I _ q2 Expressions JHU and (|100|l follow from 

^(-|,..), 

P - (0,ki), N = (0,k2), 

since k2 — k'2 = q — q' = k'l — ki = (see Fig. |5land Ea. l|^ '). Introducing the four-vector L\ = ewpaM'^ K^A'^ , 
where, 

M = N + P = (0,N + P), 

K = N-P = (0,N-P), (102) 
A = k[-ki = {Vt,0), 

Q can be represented as 

Q-M Q K Q ■ A , Q ■ L , , 

^ M+^K+^A + ^L. (103) 



M2 • ' A2 ' i2 

Eq. HH7|I can then be decomposed in the invariants Pi of Eq. (|89|l as |2 



J2[pf{w,t) T (-l)Vf (4to2 -w- t,t)]P, = ^^L^ 

i 

X [{A^ + 2mC^ + to2£±) Pi + (6± + mV^) P2 + f^Pg - T^Pi] , (104) 



A Nucleon Pole Contribution for Pseudo-scalar Coupling 
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with 



M2 



^2 



d^QB, A, ^-j^^ — 

(Q • M)2 



1 



1 



± _ 



d^Q B^'^ B^ 



{2K^ + A/2) 



M2 

(g ■ M)2 

M2 

(Q.M)2 ^ (Q.i^)2 



2i^2 - Af2) 



(Q-K) 



X2 



(105) 



Here Af2 = uu and — t — Am? — w — t. Note that the Bjorken-DreU metric used here leads to some sign differences 
with respect to references ^3 1^ where the Pauh metric was used. 



A. Nucleon Pole Contribution for Pseudo-scalar Coupling 



If the invariant amphtudes A{B)^^^ are expressed as fixed-t dispersion relations j52| | one can write the pf^ given by 
Eqs. pU4|l and pU5|l in the spectral form Following Ref. [s^l one can use the Cutkosky rules [s^ to calculate 

the imaginary part of the pf from Eq. pU5|) . The expressions for the yf{w^ t) of Eq. H94|l in terms of double integrals 
over si and S2 on quantities depending on elastic pion-nuclcon scattering absorptive parts, are given by Eqs. (2.6) to 
(2.8) of Ref. ,25]. The double spectral functions yfpj{w,t) of the box diagram for pseudo-scalar (PS) pion-nucleon 

coupling can then be obtained from these equations. In this case the A{B)^^^ pg for the ttN nucleon Born term are 
(see e.g. Eq. (A.6.30) of Ref. |52]), 



/1± 

^1(2), PS 



-^1(2), PS -9 ( ™2 _ 



1 



1 



"T- - Si(2) - *e "T- - si(2) - *e 



(106) 



The absorptive parts of ^f(2) PS 'L^ro and those of ^^2) PS given by iig (5(si(2) ^ w ) for the direct ter 
leading to, 



7r2 / g-^X^ K{w,t) 



±- 



2 \4ttJ (4^2 - i- w)2 
( g^\^ m K{w, t) 



2 V47r/ ^t{Am'^ -t-wf 



2 \Ati) Vt 

0, 



~2w -t 



4to2 _ i 

"iy2is:(u;,i)2 



1 



(i - 2m') 
+ (t-2/i2)' 



w (4r7i2 — t — w) K{w, t)2 



(107) 
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with 



K{w, t) = w-^'^ [{t ~ V) (w - W) - 4/i 



41-1/2 



We define, 



A{v,t') 



+ CXD 



dw' 
w' — y 



K{w',t'), C{y,t') 



-\-oo 



dw' 

{w' - yY 



K{w',t' 



(108) 



(109) 



where Xq is the lower limit of the w' integration such that K{w' , t') is defined, i. e., (t' — 4/i^)(w' — 4m^) — 4/i'' > 0, 
which leads to, 



With the definition. 



we have \5t 



with 



r — 4/i-^ 



r]{w',t') = 4/ - {t' - A^j.^){w' - 4m2) = (.xo - w')it' - 4^i^), 

+ C(4m2-t',i') 



nfg'^Ym^ {t'-2p'^f \ A{w,t') - AiAm"^ - t' ,t') 



Am? — w — t' 



p±(z«,i') = ±f (l^j ^ [^i^(0,t')+-^2A(«;,t')+-^3C(4m2-t',O+^4^(4m2-t',t')]. 



The Ti are given by. 



~ w (4m2 - t') 



^2 = - 



77(w,i') (2w + i' -4m2) 



4m^ — w — t' ' 



w {w + t' ~ Arn^Y 



The expression for p^iw, t') is 



with 



P2l>(«^,t') = ±1(1^) ^(t'-4m2) [giyl(0,i') + a2A(,i;,t')+e3C(4m2-i',i') + e4^(4™'-t',t')] 



and. 



Si = 

^3 = 



^1 



4m2 - t' ' 

•^3 

4to2 _ t' ' 



^2 = 



^2 



2u; + i' - 4m2 ' 

Ga — —Qi — Q2- 



(110) 

(111) 

(112) 



(113) 



(114) 



(115) 



(116) 



(117) 
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This square root divergence will disappear after performing the integration over t' . 
The integration over t' in Eq. H95() is performed analytically. Let us define, 



41-1/2 



(127) 



A{w',y) 



dt' 



K{w',t'), C{w',y) 



dt' 



K{w',t' 



(128) 



where xq is the lower limit of the t' integration such that K{w' , t') is defined, i. e., {t' — 4:fi^){w' — Aw?) — AjjL^ > 0, 
which leads to, 



4^4 



(129) 



One obtains, with y = Aw? — w — w' and w' = Aw^ — w' , 

afiw',y) = (l^) [aii(w',y) + ai2(w',y)], 
^ ^ ^ '"aii(w',y) 



In Eq. fnn|) . 

aii{w',y) 

ai2{w',y) 
ai2{w',y) 



a,iw',y) = ±-^- 
a^{w',y) 



+ ai2{w',y) 



(130) 



TO"' 



2 \ 2 -, 



04 = ±77 TT 



2 \47r/ wy/v? 



ri{w' , y)A{w' , y) — rj{w' , w')A{'w' , w') 



{y - 2^1^)^ A{w', w') - {w' - 2^1^) {y - w ~ 2(1^) A{w\ w') + w {w' - 2fi^yC{w\ w') 

{w — w')ri{w' ,y)A{w' , y) — [w'^ + {w' — w)r]{w' ,y)]A{w' , w') + ww'ri{w' , w')C{w' ,w') 
{w + w')rj{w' ,y)A{'w' , y) + [u? — (tu + w')ri{'w' ^y)\A{w' , w') + ww'r]{w' , w')C{w' ,w') 



The functions A{w' ,y) and C{w' ,y) with y < are. 



A{w',y) = 



In- 



\/-viw',y)y 1 _ ./Ze 
V V-- 



(131) 



(132) 



C{w',y) 



w' {y - 2m'^) + 2ii^ 



yq{w',y) yr]{w',y) 
We have calculated the discontinuities of the pf along w, i. e.. 



A{w',y). 



(133) 



1 f^°° dt' 

lim [pi{w + te,t') -pi{w - ie,t')] = 2ilmu, pi{w,t') = 2i— / " 

e^O 77 ^4^2 t - t 



Im^, pf{w,t'). 



(134) 
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We introduce the following explicit notations 

= Rept{w,t'), 

pf"°{w,t') = T{-iypf''°''{w^^m^~t'-w,t'), (135) 

Ini„ pf '^™(u;,i') = Im„ pf{w,t') = y^j^{w,t'). 

Here again, the last equality of Eq. p35|) follows from the application of Eq. H9t)|) . In the expressions pi3|l . (|116|l . pi9f) 
and (|124|l . the imaginary part of the pf{w,t') comes only from A. There are only three A terms: A{'w,t'), A{0,t') 
and A(4m^ — t',t') in Eqs. H113(l and pi6|l that have a non-zero imaginary part for y > xq. However, from Hll()|) 
one cannot have 4m^ — t' > xq, so, a fortiori one cannot have 4m^ — t' — w > xq and only the box diagram has 
an imaginary part which justifies the superscript "box" in the expression for the imaginary part. The superscript 
"cro" denotes the contribution from the crossed diagram. We write 

Hf' 

pf^°^w,t) = 3 / -— pf^o^w,t'), (136) 

fit' 

pf"°{lt) = 3 / -_p±-(f,t') 

/ — ^a^iw^Am^ -w-w'), (137) 

dt' 

lm^pf{w,t) ^ ^ I —-yf^^{w,t'). (138) 

It can be seen from studying the high t' behavior of the functions (|135|l that they decrease like 1/t' or \\i{t')/t' for 
t' — > +00, so that the integrals (|136|l . (|137|l and (|138|l on t' converge. A similar convergence holds for the integrals 
on w' . Note that the physical values of t{t) being negative, there are no poles in the t'{w') integration, so there is 
no further imaginary part. In the present form, the numerical calculation of Eqs. H13t)|) . H137|) and (|138|l requires one 
numerical integration over t' or w' . 

We have seen (Eq. H126|l ') that A(ij,t') has a divergence in y ^ x^ . In the expressions for integration on t' , y takes 
the values 4m^ — t' , 0, 4m^ — t' — w and w (see Eqs. H113|) . (I116|l and H135|l ). The divergence at y — xq is in the 
expression for p^^ in Eq. (|113|l where the first argument of A is w, so that only the box diagram has this divergence. 
Fixing w and using Eq. H11(J|I the condition w — xq gives for pf^°^{w, t') a discontinuity at t' — tg with 

*o = V + ^^. (139) 

The integral 1)136(1 can be split into two pieces, 

pf^^^w.t) = i / pp™(«;,0 + i / pf ^°'^(^,0. (140) 





r + 00 




'4^2 t' 




p + 00 




'4^2 t 


T(- 


TT 




p + CC 




'4^2 t' 



TT ,/4^2 t'-f TTJt, t'-t 



The first integral, of the type. 



*° dt' 



4/^2 y/t'o - t' 



F(t'), (141) 



becomes, with x = — f'. 



2 / F(t^ - x^). (142) 
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This transformation is applied to F{t') = \/t'Q — t' / {7r{t' ~ t)) pf^°^{w,t'). The second integral can be recast into 
a finite domain by the following change of variable 

, \{x + 1) - 2t'Q 



X - 1 



where A is a free parameter. One then has to calculate, 



(143) 



For the calculation oi pf"°{t^t) and Im^ pf{w,t) there is no discontinuity in Iq. For pf"°{t,t), one can transform 
the first term of (|137|l into an integral of the type H143() . For the imaginary part the domain of definition of K{w, t') 
restricts the integration interval to [Iq, +oo]. There is a singularity at the lower limit in Xj ^t' — t'^ which we treat as 
above (see Eqs. (|141|) and (|142|l 'l. 

We have (see Eqs. (fTU7| . (fn5|) and (HSU)), 



The proton- proton amplitude being a pure isospin one state (1=1), one has (see Eq. HSyfl). 



(144) 



M'-^ = {V't^°"{w, t) + Im,p+{w, t) + 5p+"°(t, t)] P,. 



(145) 



1=1 



We give below the formulae which express, for a given isospin state, the nucleon-nucleon helicities amplitudes Lpi in 
terms of the pi amplitudes. One has 



= 

Vi = 

ip4 = 



-— (pi - 2mDp2 + m^D^ps) + [D — 



Pi, 



1-z f E 



2E^ 



2 \TO 

l + z 



-P2 + E-p3 + P4 H , 



{pi - 2mDp2 + m^D^ps + Dp^ 



(146) 



1- z fE^ 2E'^ 



^fi P2 + E-p3+p4 ^ 

TO^ TO TO^ 



P5 



EVT^ 
2to 



2^' 2 

-Pi H P2 - m Dps - _P4 , 

m 



where D = {E"^ + k'^)/m? and z = cos 6* (let us remember that in the center of mass of the nucleon-nucleon w-channel, 
w — AE"^ = 4:{k'^ + TO^) and t ~ — 2fc^(l — z)). The Saclay amplitudes a, b, c, d and e are given in terms of the (pi 
amplitudes by the following expressions [3^ . 



(Vl + tp2 + f3 - '^i) z ~ 4(^5 v^l - 



- (.^1 - ip2 + ip3 + <P4), 
1 / 

-(-(pi + ip2 + (P3 + 'PAj, 
i {ifl +(p2-f3 + (fid, 



(147) 



{(fil + ip2 + ^3 - fi) \/\~ z"^ + 4<^5Z 
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B. Nucleon pole contribution for Pseudo-vector coupling 

For pseudo-vector pion-nucleon coupling, the yl(_B)^2) py ^'^^ ^^'^ nucleon Born term are (see e.g. Eqs. (A. 8.1) 
and (A.8.2) of Ref. $^), 



4+ — 4- — n 

^1(2), py ~ ^i(2),py ~ 

D+ 0+ r?— r>— ^ 

^1(2), PV ~ ^1(2), P5' -^1(2),PV ~ -°1(2),P5 2m^ ' 



(148) 



One can then use Eqs. (|104|l and pU5|) to calculate the corresponding pf. One recovers first the contribution for pseudo- 
scalar coupling from B^pg B^ pg given in the above section. There are further bubble [bub) diagram contributions 

arising from A'^ py py and proportional to ( — j and from B^ py B^ py and proportional to ( ^ — - j . There 

2 

are also triangle [tri) diagram contributions arising from B^ py B^ py and proportional to ^ ^ ^1(2) PS' finds 



2m2 1(2). 



327rm2 V t 



.4 I f _ 4^2 ^> _ 4^2 

Wtt V - 2^2 



5^ t^-V (149) 



2567rm2 



P4-*"(i') = tJ^a '-^^ (*'-V) [/o(i')-/2(0]- 



In the above formulae (|149|) 



/2„(0= 1^'^^^ h{t') = ^[-2 + aIo{t% (150) 



with fi^ - t'/2. For < 4m2, p = if] = i2[{jn^ - t'/4){t' /A - p'^)Y''^ and 



2 /3 

/o(t') = -arctan-. (151) 



For t' > 4m2, (3 = 13 = 2[{t' /A - m^){t' jA - Ai^)]i/2 and 

4(0 = ^ In 



a + /3 



a- ^ 

The functions /o(i') and /2(i') are continuous at t' = Am} and 



Furthermore P2 '^'(i') is also continuous at t' = Am , with, 



(152) 



/o(4m2) = -, hiAm^) = ^. (153) 
a 6a 



W)^3W) ^ 2im^-p^) 
t'_4m2 t' - 4m2 (^2 _ 2m2)3 ^ ^ 
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For t' oo the functions /o,2(i') — > (2/^') \n{t'/m^) and (t'/4 — ^^)[/o(t')~^2(^')] ^ constant. It can then be seen from 
Eqs. (ESIl and (|149|l that the amphtudes p^*"{t) and P2*"{t) are convergent, Pi'"^''{t) and p1*"{t) are logarithmically 
divergent and p^^^''^[t) is linearly divergent. Introducing an upper limit of integration, t^/, in the integrals (|93(l . one 
finds 



with 



-bub 9^ 

Pi " (i,tM) = -i^ofi 2 4 ^4(i, ^m), 

/i(t,iM) = /i(0,iM) + Ji(i,tAf), 

/l(0,iAf) = -2X(tM) - hn{X{tM), 1), 



Jl(i,iM) = -t[(2/3)x3(tM) -/i(0,tA/)], 



(155) 
(156) 

(157) 
(158) 
(159) 

(160) 
(161) 
(162) 
(163) 



Ji{tMi) = 2X{t„^) + X{t)hn{X{tM).X{t)), 
X(t) = UY, Z) = In , 

Ii{t,tM) = hiOjtAl) + Jl{titM) + Ji{t,tM), 
14(0, Im) — tMX^{tM) — Q{lJ''^ /tM)Jl{tM,tM), 



Jlit.tM) = t [{2/3)X''itM) + X\t)Jiit,tM)] ■ (164) 

The logarithmic divergence of Pi^^^{t) is seen in Eq. (|159(l and the linear divergence of ''"''(<) in Eq. rrn^ . 

The p^^'^^(t'), p^*"(t') contributions and those of P2^"{t'), (t') can be interpreted as coming from the 

correlated 27r-exchange in the S-wave (cr or /o(600) exchange) and in the P-wave {p exchange), respectively (see 
Eqs. (2.29), (2.30), (2.31) of Ref. ^ and Eqs. (2.12), (2.13) of Ref. However, here there is no contribution to 

p^ as A~ = 0. 

It is interesting to examine the box- and crossed-diagram contributions to the bubble and triangle diagrams arising 
from the pseudo-vector coupling case. Performing a similar reduction to that of Ppy (see Eq. (26) in section II. B) one 
can write the direct and crossed, TffirN amplitudes as 

TTf = rlrllAl^ + 7^ • Q(S., - S.)], (165) 



= rlTl\A\^ - 7^ • Q(S.-, - B,.)\. (166) 
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Here the Pauli matrices Ta(p) represent the isospin coupUng of the pion field T^aip) to the nucleon 1. In Eqs. (|165|l and 



2 ? 9 

9 n .9 „ .9 



Note that 



A^^, Br = j-, B^ = ^—. (167) 

2m Am m^ — s 



J2 ^a^/^^a^l = ('^"'3 + *ea/37'^7) ('^"/^ + ^'^o.P'y'^j) = 3 " • T^. (168) 

a, (3=1,2,3 a,f3 

The 27r-exchange box diagram contribution arises from 

= (3 - 2t1 • t2) {2A2i1i2 + {Bs, - B,,)A{l^-r^ ■ Q + 1^72 • Q) 

+2 [B,,B,, - B,(B,, + BsJ + 7I . Q72 . Q} . (169) 
The crossed diagram contribution is, 

rp rpaf3rpa(3 . rpaPrpaP 

J^CrO — J- Si -'-S2 "T^Sl ^5-2 

= (3 + 2t1 • t2) {2^2^112 + (B,, - BsJA{l^f^ ■ Q + l^f^ ■ Q) 

-2 [Ss.Bs, - Br{Bs, + B,,) + Bl] 7I • Q72 . Q} . (170) 

In Eqs. I|169|l and (|17()|l the terms Bg^Bs^ and Bg^Bg^ correspond to the 27r-exchange nucleon pole contribution for 
pseudo-scalar coupling. For the bubble and triangle diagram contributions, use of Eqs. I|l()4|l and (|l()5|l leads to 

pt\"At') = T pt'::oit') = tI p^'-in m) 

It can be seen that the total pfl]^'^^*^"\t' ) + pfcr'o'*"\^' ) from the above formulae 1171(1 are in agreement with the 
formulae p49|l calculated from Eqs. (|148|) . Formulae (|171|l satisfy the relations of Eq. p35|l between the box and 
crossed pi. 



C. Comparison with Feynman diagram calculation 

We have checked, at different energies and angles, that the Saclay amplitudes given by Eq. (|147|) for both box and 
crossed diagrams for pseudo-scalar coupling calculated in Section II-A from the Feynman diagram expressions with 
no form factor (A 00) are in agreement with the corresponding dispersion relation results Eq. p36|) . Eq. p37|) 
and Eq. (|138|l . We found that the short range contribution of the crossed diagram coming from the term containing 
af (i(;',4m^ — w — w') in Eq. I|93() is crucial for the agreement. Numerical comparisons will be shown in the next 
section. 
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For the pseudo-vector coupling, the comparison of the bubble and triangle contributions Eq. H149|) is less straight- 
forward due to the divergence of some of these terms when integrated on t' in the dispersion relation or of the 
corresponding terms in the limit of A ^ oo in the Feynman method. The results of both approaches for the conver- 
gent invariant amplitudes pf'^"{t) and P2^"{t) (see Eqs. (|149|) ') compare well. The logarithmically divergent p^^^^{t) 
and p'^*^"{t) require one subtraction and a comparison can be made with the dp^'"''' (t) / dt and dpl*"{t)/dt. For the 
linearly divergent p^^^^{t) the second derivative, d"^ p^^"^^ (t) / dt'^ , converges and can be compared with the Feynman 
calculation. 

To obtain these amplitudes in the Feynman calculation the amplitudes p4 and pi are extracted by solving for them 
from Eqs. 1] 146(1 and H147() thus removing the explicit dependence on the cosine of the scattering angle, z. Then the 
derivatives with respect to z are essentially the same as the derivatives with respect to t in the case of the dispersion 
relations allowing a direct numerical comparison. Since this process requires numerical derivatives and a calculation 
for large A it is more difficult, but a good agreement is obtained. 
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FIG. 6: Contributions to the real (left) and imaginary (right) parts of the amplitudes h, d and e as a function of energy for 
pseudo-scalar coupling. The dashed line corresponds to the box, the dashed-dot line to the crossed-pion diagram and the solid 
line to their sum. 



VI. RESULTS AND DISCUSSION 



Test calculations were made with and without the antiproton poles (po = —E). At low values of Piah or for large 
A they give significant contributions and are essential to get agreement with the dispersion relation results. However, 
for realistic values of A at higher energies they become unimportant and can be neglected. 

Figure El shows the behavior of the real and imaginary parts of the amplitudes as a function of incident momentum 
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FIG. 7: Comparison of the real (left) and imaginary (right) parts of the amplitudes b, d and e corresponding to the pseudo-scalar 
and pseudo-vector couplings. The dashed curve represents the PS and the solid curve the PV. 



for the pseudo-scalar coupling. For the crossed-pion diagram only terms vi^ith a single nucleon or pion on shell and 
the principal value integral on survive. That is to say, only the contributions from the go integral need to be 
treated specially, the rest of the indicated integrals have no singularities. For the box diagram, in addition to the two 
contributions just mentioned, the diagram in which two nucleons are on shell also gives a real contribution (in the form 
we are looking at now, it becomes an imaginary contribution to the amplitude after application of the factor of i to 
be consistent with the dispersion relation approach). That is, the integral on the magnitude of the three-momentum 
also has a (5-function part which contributes. 

In figure ini (imaginary part) we see that the principal- value part of the integral goes rapidly to zero for the crossed 
diagram as the momentum drops below 500 MeV/c. This behavior can be traced to the form factor and its dependence 
on go- For low energy the form factor becomes independent of go? since the pole moves far from the real axis. Since 
it is the only pole off of the real axis the integrand limits to a function with only poles on the real axis and hence 
contains only the (5- function parts. 

The imaginary part of the box contains both the principal value and double pole contribution. At low energies the 
real parts of the amplitudes are dominated by the term with two nucleons on shell. By 2 GeV/c that situation has 
been reversed with the crossed-pion principal value dominating. 

Figure [3 compares the results for the pseudo-scalar and pseudo- vector couplings. One sees that the real part of the 
Saclay amplitude b is very large for the PS coupling but is greatly reduced for PV coupling. The d amplitude shows 
only a modest difference between PS and PV coupling while the e amplitude is again significantly modified. We note 
that the h amplitude is a mixture of central and spin-spin character, the d amplitude is of tensor character and the e 
amplitude is of spin-orbit character [3| ■ 

Figures IHl and 1^1 show a comparison of the box and crossed-pion contributions. It is readily seen that the crossed- 
pion diagram dominates at high energy. This result can be traced back to the relative ability of the two diagrams to 
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FIG. 8: Results for the crossed (left) and box diagram (right) for the b amplitude. The dash-dot curves show the results of the 
dispersion relation calculation, equivalent to the Feynman calculation for A — > oo. 
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FIG. 9: Results for the crossed (left) and box diagram (right) for the d and e amplitudes. The curves have the same meaning 
as in Fig. |S| 



transfer energy and momentum. Because factors of momentum transfer in the numerator cancel the decrease coming 
from the propagators, the integrand without form factors does not fall off with increasing momentum transfer so 
that the decrease in cross section is mainly due to the form factors. Since the form factor falls off with increasing 
pion momentum in the center of mass of the nucleon, one can see that for the box diagram with the four-vector q 
connecting vertices with {E, k) and [E, — k) both momenta in the respective centers of momentum cannot be zero (or 
even small) at the same time. The crossed diagram does not suffer from this contradiction since q connects — k and 
k' and, for some value of the 4-vector, q, there is a chance to minimize the momenta at both vertices. We now look 
at the behavior of the form factor chosen for these calculations. 

The crucial element for the understanding of the behavior of the box and crossed diagrams is the argument in 
Eq. Since one nucleon is always on-shell we can write as 



(172) 



where a = ±k or ± k' and b — q or q' . 

The order of magnitude oiq~ |q| will be the same as |k| since g cannot be taken as small to maximize the integrand 
because of the powers of q in the numerator. For values of momenta greater than the mass of the nucleon, the first 
term, (a • 6)^, will dominate so we need only consider the size of that term to obtain some sense of the behavior of 
the form factor. Out-of-plane values of q (in the y direction in the system we are using) only increase Z^ so consider 
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Qy = 0. Let X be the cosine between the incident direction k and the vector q and consider the pole at qo = lu in an 

extreme relativistic limit E ^ k {k = and uj q. 

It is useful to look for the minimum in the four values of in order to maximize the form factor. Corresponding 
to the four factors in Eq. (|85|1 we have (neglecting the second term in Z^) the following four values for a • & for the 
box diagram; 
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1) Euj~k-q-y kq{l - x); 2) Eto + k-q^ kq{l + x); 



3) Eu-k! ■€{ ^ k[q{\ - \/\ -a;2) + fc]; 4) Eu + k' ■({ ^ k[q{\ + - x"^) - k]. 



(173) 



Since the first and second vertices involve the sum and difference of the same quantities, they cannot both be small 
at the same time. The form factor being an even function of x, there is an extremum at a; = and, in fact, it is a 
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maximum. If we choose x = and the fourth value of a • 5 to be zero with q — then, in this limit 

1 1 



•T-max 
^hox 



1 



(174) 



We can apply the same considerations for the crossed-pion graph from Eq H86(l . The values of a • 6 at the four vertices 
are 

1) £;w - k • q' ^ kq{l - x) - 2) + k • q ^ kq{l + x); 



3) £;w - k' • q ^ fcq(l - ^/l-x'^)] 4) + k' • q' ^ fcg(l + ^Jl-x^) - k^ 



(175) 



It is possible to find values of q and x such that any of the values of is zeroed, unlike the case of the box diagram 
where the third value can never be zero. For example, one could take a; = — 1 and q = ^k to set the first and second 
values of to zero. The maximum comes about, however, if we take the first and fourth values of a • 6 to be zero 
which leads to 



g=(2-\/2)fc; x = -l/V2. 



(176) 
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The second and third values of a • 6 become 

(3 - 2%/2)P w 0.1716P, (177) 
which is much smaller than what one was able to achieve in the case of the box diagram. The form factor is then 

(178) 



'^crossed 



0.1716fc^ 
mA 



Thus, not only is it possible to make two of the arguments zero at the same time but the other two are relatively 
small also. Numerical studies of the relative sizes of the form factors alone confirm that the values of the maximum 
of the crossed form factor occurs at the values given by Eq. H176|) . While it is necessary to go to very high energy to 
justify the ultra-relativistic limit used above for illustration, at PLQb=5.5 GeV/c typical values of the crossed diagram 
form factor are more than two orders of magnitude larger than that for the box diagram. 

Since the crossed-pion form factor peaks at definite values of q and x at high energies, analytical predictions of the 
relative size of the different amplitudes can be made. From these values we find for the Saclay amplitudes, in the very 
high-energy limit 



-2h- d^O, 



(179) 



which was verified by a calculation at P^afc — 46 GeV/c. These values lead to a value of Cnn which is again 1/3. In 
the energy range where we compare with data we are far from this limit, however. 
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We see that the values of the spin separated cross sections shown in Fig. ^]indicate a preference for the PV coupling 
and values of A in the range 1.2 to 1.4. Since A has been limited from other sources (see Section liVA|l . the agreement 
is obtained in a nearly parameter free manner. 

Some experimental values of the Saclay amplitudes are known at high enough energy that unitarity corrections 
can be expected to be small enough that they can be compared with the present calculation. Arndt et al. jH3 give 
amplitudes up to Piab = 3.82 GeV/c and Bystricky et al. |53| present several single energy values. Ghahramany and 
Forozani found values of the amplitudes at two energies in reasonable agreement with Refs. |53| and js^- Figure 
1121 shows a comparison of the present calculation with the SAID values (broken line) and the Bystricky et al. points. 
The agreement for the pseudo-vector coupling is seen to be considerably better than the pseudo-scalar one and gives 
a reasonable representation of the data except for the imaginary part of b. 
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One may ask if the difference in fall-ofF of the crossed and box diagrams is a general result of the kinematics or if 
it depends on the particular form of the variable used in the form factor. To attempt to give a partial answer to 
this question we have calculated the pseudo-scalar result for the crossed diagram with a different, though superficially 
similar, variable. Haberzettl et al. [43 use a variable which is of the form 



W^^(Pl,P2,P3) 



Am? 



(180) 



where we have chosen the normalization such that the limit for large values of the invariant masses the limit matches 
that of Eq. H8UI) . Figure E| shows the comparison of the two calculations and one sees that the fall-off for this variable 
is very similar to the presently used form factor. The rise at small momenta is at least partially due to the contribution 
of the principal value integral since the pole in the complex plane does not move far from the real axis as the energy 
goes to zero so that causality is not respected in that limit. 

It is useful to examine the decrease of the cross section with increasing energy. The prediction of slope of Brodsky 
and Farrar (see Fig. ^ based on counting of internal propagators and our calculation of 2|5p are very similar. 
The high energy limit of a form factor has been compared with the picture of the propagators between interacting 
constituents by the authors of Ref. |^ for the case of scattering from nuclear constituents. The two-pion-exchange 
diagram may be able to provide a bridge between the low and high-energy points of view. 

The present result bears directly on the question of the energy regime where the transition from a color singlet 
hadronic exchange to a quark-gluon basis might reasonably take place. This, in turn, may impact the question of color 
transparency (see, for example Jain et al. [BJ). If the exchange of color singlets (pions in this case) continues to be 
important through 12 GeV/c it might negate the basis for color transparency (dominance of quark-gluon exchange) in 
the moderate energy range P^afc ^ 12 GeV/c. This is the entire range covered in the color transparency experiments 
of Carroll et al. However, it may be that, even if color singlet exchange continues to be very important for elastic 
scattering, quark-gluon exchange might dominate the inelastic processes. 
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The comparison with the spin transfer observable Cnn is given in Fig. ^]for several values of A. For the pseudo- 
vector choice of coupling it is seen that the qualitative features are reproduced. One, in fact, needs a mechanism (or 
several mechanisms) with considerable structure in order to reproduce the data. The plateau in the data from 4 to 8 
GeV/c is not found but there is a minimum in the calculation in this region. 
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VII. SUMMARY AND CONCLUSIONS 

We have calculated the contribution of the box and crossed two-pion-exchangc diagrams to proton-proton scattering 
at 90°^ for laboratory momenta up to 12 GeV/c. The cases of both pure pseudo-scalar and pure pseudo-vector 
pion-nucleon coupling are treated. 

We assume that the interaction of the pion is with the valence quarks within the nucleons and use an effective 
field theory obeying Feynman rules to describe the composite system. At each pion-nucleon-interaction vertex we 
introduce a relativistic Lorentz-invariant form factor. These form factors are related to the convolution of the nucleon 
and pion sizes and represent the pion source distribution based on the quark structure of the hadronic cores and 
explicitly introduce an interaction range of this size. These form factors are functions of the four-momenta of the 
exchanged pions and scattering nucleons. One can use any two of the three momenta which converge at a given vertex 
to calculate the function so that the dependence on the nucleon and pion momenta displays a high degree of symmetry. 
This behavior can be traced to the fact that the basic scalar on which the form factors depend can be written as the 
square of the four-dimensional cross product of any two of the three four- vectors. 

While we believe the form factors that we introduce are quite reasonable, we make no claim as to their uniqueness. 
However, there exist a set of conditions, based on short reaction time, causality etc., that any relativistic Lorentz 
invariant form factor should satisfy. Our form factors obey these rules, but other form factors used by some authors 
do not. 

As a check of the calculation, we compare the two-pion-exchange Feynman-diagram amplitudes for both pseudo- 
scalar and pseudo- vector coupling calculated with point-like nucleons with those obtained from the dispersion relation 
approach and find agreement. The antiproton poles in the Feynman calculations are essential for the comparison 
but for a more realistic calculation with a finite value of the form-factor range. A, they are unimportant. While 
performing this study, we found it essential to include a contribution of the crossed diagram which was neglected in 
the dispersion-relation calculations of the two-pion-exchange Paris potential because of its very short range. 

The numerical technique for calculating the 4-dimensional Feynman integrals, taking into account multiple poles is 
presented. Standard transformation methods are not always applicable in the presence of form factors. 

In the Feynman diagram calculations, using the form factors constrained by the valence quark distribution, com- 
parison with experimental data favors the pseudo-vector coupling over the pseudo-scalar one. This is seen in the 
comparison with the magnitude of the spin-separated cross sections, in the agreement with amplitudes extracted up 
to 4 GeV/c and in the spin transfer parameter Cnn, although in this later case it is only the general behavior which 
is correctly given. 

The exchange of one pion is important only at momentum less than Piaf,=l GeV/c. While the strengths of the box 
and crossed diagrams for the exchange of two pions are comparable for laboratory momenta below 2 GeV/c, for larger 
momenta the crossed diagram dominates, mainly due to the mathematical structure of the form factors and to the 
fact that the kinematics of the crossed diagram allows a repartition of momenta at the vertices in a favorable manner. 
An important contribution arises from the principal-value contribution of the integrals which is non-zero when form 
factors are included. 

While here we compare only with the extracted spin-separated cross sections at 90° ^ , our calculation can be 
extended to all scattering angles. Further studies such as unitarity corrections can be considered. 

To conclude, we have shown that the two-pion exchange plays a significant role in proton-proton scattering even 
up to Pi,ab=12 GeV/c. This result suggests that the importance of the exchange of color singlets in elastic scattering 
may extend to higher energies than expected. 
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APPENDIX A: PROJECTION ONTO THE SACLAY AMPLITUDES 

The Saclay amplitudes [s^ are defined by the equation 

M{kf,ki) = ^[(a + &) + (a - b)(Ti ■ n(T2 ■ n + (c + d)(Ti ■ m(T2 • m, 

+ (c- c?)(Ti • 10-2 • 1 + e(cri +0-2) • n] (Al) 
where the center-of-mass basis vectors are: 

k/ + ki k/ — kj ki X k/ 

^^WTT^-y "^^^^7^' "^R^- ^^^^ 

Observables at 90° are 

a = i(2|6p + \d\' + lep); aC^^ = i(-2|6p + |dp + jep); (A3) 

ai^ATjv = + lep) = (t/^^at; aCMAfMAf = ^(2|6p + ~ \e\^). (A4) 

The expressions for the amphtudes in terms of the singlet-triplet matrix elements are: 

a = i(Mii +M00 -Mi_i); & = J(Mii + M,, + Mi_i); (A5) 



c= i(Mn-M,, + Mi_i); d = yJ—{Aho+Moi); e = ^(Mio - A/qi). (A6) 

^ V2sm0 V2 

The two-pion exchange box and crossed amplitudes can be written in the form 

Af(k/, k,) =« (G^ + iai ■ H1)(G2 + zcr2 • H^) » (A7) 

where the double brackets indicate integrations over the propagators and the G and H here are related to the corre- 
sponding quantities in Section II by the factor (G^ = j^G, G^ = j^G' , etc.). With our choice of coordinate 
system G and Hy are even in qy while and Hz are odd. Any totally odd quantity will integrate to zero and those 
terms have been dropped in the following expressions. In the remainder of this appendix it is to be understood that 
each term bilinear in G and/or H is surrounded by double brackets. One has 

Afi^i = G^G^ - HlHl- Ms,s = G^G^ + • H^; Mi^^i = -HIhI + iJ^iJ^ ; (A8) 
Mo.o = G^G^ + HlHl - HlHl - HIH^; (A9) 



A'io,i — ~7=[^-^z-^x ^ ^x^z ^ ^^^y ^ G^Hy]] (AlO) 
v2 



40 



Mi,o = [-HlHl - HlHl + G'Hl + G'h'J . (All) 

Prom these expressions we can construct the ampUtudes as 

a = G^G^ -HlH^; b = G^G^ + H^H^; c = -H^H^ - HIH^^; (A12) 

d=^Q{HlHl + HlHl)- e = i{G'Hl + G^Hl). (A13) 

With the symmetry due to the identical particles 

a{9) ^ a{e) - a{-K - 9), (A14) 

6(61) ^ 6(61) - c(7r - 9) = G^ {9)G^ (9) + (tt - 6l) • (tt - 6*), (A15) 

c{9) = -b{TT - 9), (A16) 



"^^^^ " ^[^^'^^^^'^^^ + Hlml{9)Hl{T: - 9)Hl{n - 9) + HUtt - 9)Hl{Tr - 9)], (A17) 



e{9) ^ G\e)H^{e) + G^{e)H^{9) + G\tt - 9)H^{n -9) + G^{7r - 9)Hl{n - 6). (A18) 



At 90° 



a = 0; 6 = G1G2 + h1-H2 = -c; d = 2{HIhI + HIhI)- ■ e = 2{G^H^ + G^H^, (A19) 
where the quantities are evaluated at 90°. 

APPENDIX B: CALCULATION TECHNIQUE 

Consider the interior integrals over q and go with the integrals over the solid angle to be done at the exterior. We 
apply two formulae depending on the case. In the integral of the type 

there are two possibilities. If xi and X2 are distinct for all values of y then we can write the integral as 
J J [x - xx{y)][x - X2{y)] J J [x - x-i_{y)][x - x^iy)] 



where 



g{x,y) = -"^^y] f{x,{y),y) + ^^^l^J{x2{y),y). (B3) 
x-i_{y) - X2{y) X2{y)-xi{y) 
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The combination f{x,y) — g{x,y) vanishes at each of the poles so that there is no singularity in the first term. For 
the second integral we have 



g{x, y)dx 



1 



.fix2iy),y) 



-xo - xi{y)][x - X2{y)] xi{y) - X2{y) x - xi{y) X2{y) - xi{y) J x - X2{y) 



1 



(B4) 



f{xi{y),y) 



±iTT + In 



xi{y) - X2{y) _ 
For the general case of n distinct poles 

dxf{x,y) 



xq - xi{y) 
Xo + xi{y) 



fix2{y)..y) 
X2{y) - xi{y) 



dy 



n"^i[x-x,(2/)] 



= / dy 



dx[f{x,y) - g{x,y)] 



zbiTT + In 



dy 



where 



g{x,y) = }^ 



1 UjM^^iy) - xjiy)] 



f{xi{.y),y) 



and the second integral is given by 



E 



fixz{y),y) 



\ \[i^i\^i{y) - ^Ay)\ 



iin + In 



xq - Xi{y) 
Xo + Xi(y) 



Xq - X2{y) 
xo + X2{y) 

dxg{x,y) 



If two poles are not distinct (for some value of y) then this method does not work. Instead we may write 



(B5) 



(B6) 



(B7) 



(B8) 



dy 



dxf{x,y) 



[x - xi{y)][x - X2iy)] 



= / dy 



1 



xi{y) - X2{y) 



dxf{x,y) 



1 



1 



x-xi{y) x-X2{y) 



(B9) 



Each of these integrals can be done separately and a second singularity has been pushed into the y integral. In the 
cases where the pole occurs in quadratic expressions it may be more efficient to take that into account as 



dy I dx- 



f{x,y) 



x"^ - xl{y)][x'^ - xliy)] 



xliy) 



dxf{x,y) 



1 



1 



x'^ ~ x\{y) x'^ — x\{y) 



(BIO) 



APPENDIX C: INTERPRETATION OF Z'' 



One can be led to the form of Eq. (|76|l from the condition given in Eq. H78|l . Since k' = k + q (for example) 
we observe the analogous condition in three dimensions is indicative of the vector cross product and are thus led to 
consider the four-dimensional version of the cross product. If we define this cross product by the use of a totally 
antisymmetric 4-component tensor, eijke, in analogy with the 3-dimensional case, for two vectors (say a and b) the 
result is a tensor 



Tij = ^ ei]kiakbe. (CI) 

0,1,2,3 

Since this tensor has 6 independent components, it cannot be expressed as an ordinary 4- vector. It is useful to separate 
the components into two classes: one in which the zero index is free and one in which it is contained in the sum. 



Toj = ^ eojktakbi = [a x b]^; j = 1, 2, 3, 

1,2,3 



(C2) 
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Tij =aobk-boak = [aoh-boa\k, i,j,k = 1,2,3 and cyclic. (C3) 
Contracting this tensor with itself with the standard metric tensor Qij = giSij; go = 1; gi = —1, i = 1, 2, 3 we find 

I Yl 9ii'9j3'TijTi,j' = IY1 9ii9jjTijTij = (aob - boaf - (ax b)^ = (a • bf - a^b^. (C4) 

The last identity may be verified by direct evaluation and has the form used for showing that it corresponds to the 
contraction of a four-dimensional cross product. 



[1] A. Lin et al., Phys. Lett. 74B, 273 (1978): Energy dependence of spin-spin forces in OOlrn elastic pp scattering. 

[2] T. S. Bhatia ct aL, Phys. Rev. Lett. 49, 1135 (1982): Spin Correlation for pp Elastic scattering at Qc.m. ~ 7r/2 tn the 

Energy Region of Dibaryon Resonances. 
[3] E. A. Crosbie et al., Phys. Rev. D 23, 600 (1981): Energy dependence of spin-spin effects in pp elastic scattering at 

9(f cm.. 

[4] G. R. Court et aL, Phys. Rev. Lett. 57, 507 (1986): Energy Dependence of Spin Effects in p| +p| — > p + p. 

[5] D. G. Crabb et al, Phys. Rev. Lett. 41, 1257 (1978): Spin dependence of High-p\ Elastic p-p Scattering. 

[6] J. R. O'Fallon et al., Phys. Rev. Lett. 39, 733 (1977): Spin-Spm Interaction in High-p\ Elastic p-p Scattering. 

[7] C. W. Akerlof, R. H. Hieber, A. D. Krisch, K. W. Edwards, L. G. Ratner and K. Ruddick, Phys. Rev. 159, 1138 (1967): 

Elastic Proton-Proton scattering at 9(f and Structure within the Proton. 
[8] R. C. Kammerud, B. B. Brabson, R. R. Crittenden, R. M. Heinz, H. A. Neal, H. W. Paik and R. Sidwell, Phys. Rev. D 4, 

1309 (1971): Large- Angle Proton-Proton Elastic Scattering at Intermediate Momenta. 
[9] K. Abe, B. A. Barnet, J. H. Goldman, A. T. Laasanen, P. H. Steinberg, G. J. Marnier, D. R. Moffett and E. F. Parker, 

Phys. Rev. D 12, 1 (1975): Wide-angle differential cross sections for elastic proton-proton scattering in the region of 

A(1236) production. 

[10] S. J. Brodsky, C. E. Carlson and H. Lipkin, Phys. Rev. D 20, 2278 (1979): Spin effects in large-transverse-momentum 
exclusive Scattering processes. 

[11] G. R. Farrar, S. Gottlieb, D. Sivers and G. H. Thomas, Phys. Rev. D 20, 202 (1979): Constituent description of NN elastic 

scattering observables at large angles. 
[12] G. P. Ramsey and D. Sivers, Phys. Rev. D 45, 79 (1992): Spin observables for NN NN at large momentum transfer. 
[13] G. P. Ramsey and D. Sivers, Phys. Rev. D 47, 93 (1993): Proton-proton elastic spin observables at large t. 
[14] T. E. O. Ericson and M. RosarClot, Nucl. Phys. A405, 497 (1983): The deuteron asymptotic D-state as a probe of the 

nucleon-nucleon force. 

[15] J. L. Friar, B. F. Gibson and G. L. Payne, Phys. Rev. C 30, 1084 (1984): One-pion exchange potential deuteron. 

[16] J. L. Ballot and M. R. Robilotta, Phys. Rev. C 45, 986 (1992): Pionic values for deuteron observables; Phys. Rev. C 45, 

990 (1992): Nonpionic effects in deuteron asymptotic observables; J. L. Ballot, A. M. Eiro and M. R. Robilotta, Phys. Rev. 

C 40, 1459 (1989): Pions in the deuteron. 
[17] W. R. Gibbs and B. Loiseau, Phys. Rev. C 50, 2742 (1994): Neutron-Proton charge exchange. 

[18] S. J. Brodsky and G. R. Farrar, Phys. Rev. D 11, 1309 (1975): Scaling laws for large-momentum-transfer processes. 
[19] M. Taketani, S. Nakamura and M. Sasaki, Prog. Theor. Phys. 6, 581 (1951): On the Method of the Theory of Nuclear 
Forces. 

[20] M. Taketani, S. Machida, and S. Ohrmma, Prog. Theor. Phys. 7, 45 (1952): The Meson Theory of Nuclear Forces. 

[21] K. A. Brueckner and K. M. Watson, Phys. Rev. 92, 1023 (1953): Nuclear Forces in Pseudoscalar Meson Theory. 

[22] M. Sugawara and S. Okubo, Phys. Rev. 117, 605 (1960): Two-Nucleon Potential from Pion Field Theory with Pseudoscalar 

Coupling; ibid 117, 611 (1960): Two-Nucleon Potential from Pion Field Theory with Pseudovector Coupling. 
[23] S. Machida, Prog. Theor. Phys. (Suppl.) 39, 91 (1967): Nonstatic Nuclear Potential. 

[24] M. H. Partovi and E. L. Lomon, Phys. Rev. D 2, 1999 (1970): Field- Theoretical Nucleon-Nucleon Potential. 

[25] W. N. Cottingham, M. Lacombe, B. Loiseau, J.-M. Richard and R. Vinh Mau, Phys. Rev. D 8, 800 (1973): Nucleon-Nucleon 

Interaction from Pion-Nucleon Phase-Shift Analysis 
[26] M. Lacombe, B. Loiseau, J.-M. Richard, R. Vinh Mau, J. Cote, P. Pires and R. de Tourreil Phys. Rev. C 21, 861 (1980): 

Parameterization of the Paris NN potential. 
[27] J.-M. Richard, Contribution a I'etude theorique et phenomenologique des forces nucleaires, These de Doctorat d'Etat Es- 

Sciences Physique, Universite P. & M. Curie, 19 June 1975. 



43 



[28] J. L. Friar and S. A. Coon, Phys. Rev. C 49, 1272 (1994): Non-adiabatic contributions to static two-pion-exchange nuclear 
potentials. 

[29] M. R. Robilotta and C. A. da Rocha, Nucl. Phys. A615, 391 (1997): Two-pion-exchange nucleon-nucleon potential model 
independent features. 

[30] M. C. M. Rentmeester, R. G. E. Timmermans, J. L. Priar and J. J. de Swart, Phys. Rev. Lett. 82, 4992 (1999): Chiral 
Two-Pion Exchange and Proton-Proton Partial-Wave Analysis. 

[31] R. P. Fcynman Phys. Rev. 76, 769 (1949): Space-Time Approach to Quantum Electrodynamics; see also Jauch and Rohrlich, 
The Theory of Photons and Electrons, Springer- Vorlag, New York 1976. 

[32] L. Y. Glozman and D. O. Riska, Phys. Repts. 268, 264 (1996): The spectrum of the nucleons and the strange hyperons 
and chiral dynamics; L. Y. Glozman, W. Plessas, K. Varga and R. F. Wagenbrunn, Nucl. Phys. A631, 469c (1998): 
Light and strange haryons in a chiral quark model with Goldstone-boson-exchange interactions; S. N. Jena, P. Panda, and 
T. C. Tripathy, Jour. Phys. G27, 227 (2001): Masses of low-lying baryons in the ground states in a power-law potential 
model with pseudoscalar meson, one-gluon and centre-of-mass corrections; C. Helminen and D. O. Riska, Phys. Rev. C 58, 
2928 (1998): n-gluon exchange interaction between constituent quarks: A. Valcarcc, H. Garcilazo and J. Vijande, Phys. 
Rev. C 72, 025206 (2005): Constituent quark model study of light- and Strang e-bary on spectra; M. Arima, K. Masutani and 
T. Sato, Prog. Thcor. Phys. Supp. 137, 169 (2000): Baryon resonances in a constituent quark model; H. R. Pang, J. L. 
Ping, F. Wang and T. Goldman, Phys. Rev. C 65, 014003 (2002): Phenomenological study of hadron interaction models. 

[33] F. Gross and Y. Surya, Phys. Rev. C 47, 703 (1993): Unitary, relativistic resonance model for ivN scattering. 

[34] P. F. A. Goudsmit, H. J. Leisi, E. Matsinos, B. L. Birbrair and A. B. Gridnev, Nucl. Phys. A575, 673 (1994): The extended 
tree-level model for the pion-nucleon interaction. 

[35] S. Kondratyuk and O. Scholten, Phys. Rev. C 59, 1070 (1999): Consistent off-shell vertex and nucleon self-energy. 

[36] J. Bystricky, F. Lehar and P. Winternitz, J. Phys. 39, 1 (1978): Formalism of Nucleon-Nucleon elastic scattering experi- 
ments. 

[37] J. D. Bjorken and S. D. Drell, Relativistic Quantum Fields McGraw-Hill, (1965), Chapter 18. 

[38] B.D. Keister and W.N. Polyzou, Advances in Nuclear Physics, Volume 20, Ed. J. W. Negele and E.W. Vogt, Plenum Press 

1991: Relativistic Harniltonian Dynamics in Nuclear and Particle Physics. 
[39] H.M. Nussenzveig, Causality and dispersion relations, Academic Press NY (1972). 

[40] C. M. Maekawa and M. R. Robilotta, Phys. Rev. C 55, 2675 (1977): Chiral symmetry: Pion-nucleon interactions in 
constituent quark models. 

[41] M. Bozoian, J. C. H. van Doremalen and H. J. Weber, Phys. Lett. 122B, 138 (1983): Are There Bound A=2 Hyperon 

States?. 

[42] S. Nozawa, B. Blankleider and T.-S. H. Lee, Nucl. Phys. A513, 459 (1990): A Dynamical Model of Pion Photoproduction 
on the Nucleon. 

[43] H. Haberzettl, C. Bennhold, T. Mart and T. Feuster, Phys. Rev. C 58, R40 (1998): Gauge-invariant Tree-level Photopro- 
duction Amplitudes with Form Factors. 
[44] M. Gourdin, Phys. Rept., 11, 29 (1974): Weak and electromagnetic form-factors of hadrons. 

[45] M. Lissia, M.-C. Chu, J. W. Ncgclc and J. M. Grandy, Nucl. Phys. A555, 272 (1993): Comparison of hadron quark 

distribution from lattice QCD and the MIT bag model. 
[46] C. Alexandrou, Ph. de Forcrand and A Tsapalis, Phys. Rev. D 68, 074504 (2003): Matter and pseudoscalar densities in 

lattice QCD. 

[47] S. A. Coon and M. D. Scadron, Phys. Rev. C 42, 2256 (1990): ttA^'A^ couplings, the nNN form factor and the Goldberger- 
Treiman discrepancy; Phys. Rev. C 23, 1150 (1981): Goldberg er-Treiman discrepancy and the momentum variation of the 
pion-nucleon form factor and pion decay constant. 

[48] M. Bolsterh and J. Parmentola, Phys. Rev. D 39, 1304 (1989): P-wave pion m the Los Alamos soliton model. 

[49] G. Ramalho, A. Arriaga and M. T. Pena, Phys. Rev. C 60, 047001 (1999): Two-pion exchange and strong form factors in 
covariant field theories. 

[50] B. D. Keister and L. S. Kisslinger, Nucl. Phys. A326, 445 (1979): A relativistic study of pp ^ d tt"*": Forward scattering; 
G. Wolf, Phys. Rev. 182, 1538 (1969): Single- and Double-Pion Production by One-Pion Exchange and a Comparison with 
Experimental Data between 1.6 and 2.6 GeV/c; H. P. Durr and H. P. Pilkuhn, Nuovo Cim. 40, 899 (1965): Kinematical 
form factors in the peripheral model; J. Beneckc and H. P. Durr, Nuovo Cim. 56, 269 (1968): A Relativistic Model for 
Kinematical Form Factors; L. Mathelitsch and H. Garcilazo, Phys. Rev. C 33, 2075 (1986): ReUitivistic one-pion exchange 
potentials. 

[51] S. Mandelstam, Phys. Rev. 112, 1344 (1958): Determination of the Pion-Nucleon Scattering Amplitude from Dispersion 

Relations and Unitarity. General Theory. 
[52] G. Hohler, tt N scattering in Landolt-Borstein, New Series, Vol. 9b (1983). 

[53] G. E. Brown, A. D. Jackson, The Nucleon-Nucleon Interaction, North-Holland/ American Elsevier (1976) chap. XIV. 
[54] K. Nishijima, Fields and particles, W. A. Benjamin Inc. 1969, pp 419-421. 

[55] R. E. Cutkosky, J. Math. Phys. 1, 429 (1960): Singularities and Discontinuities of Feynman Amplitudes. 

[56] H. Zytnicki, Parametre de correlation de spin proton-proton a 9CP jusqu'd 12 GeV/c et echange de deux pions a 



44 



I'aide des relations de dispersion, Rapport de stage de D.E.A. de Physique Quantique, Uni. P. & M. Curie, Paris, 

http: / /Ipnhe-theorie. in2p3.fr/TH-Rapp. html 
[57] SAID program at http://gwdac.phys.gwu.edu R. A. Arndt, I. I Strakovsky and R. L. Workman, Phys. Rev. C 62, 034005 

(2000): Nucleon Nucleon Elastic Scattering to 3 GeV. 
[58] J. Bystricky, C. Lechanoine-LeLuc and F. Lehar, Eur. Phys. J. C4, 607 (1998): Direct reconstruction of pp elastic scattering 

amplitudes and phase shift analyses at fixed energies from 1,80 to 2.70 GeV. 
[59] N. Ghahramany, G. Forozani, Phys. Rev. C 61, 064004 (2000): Direct reconstruction of pp elastic scattering amplitudes at 

1.8 and 2.1 GeV. 

[60] R. D. Amado and R. M. Woloshyn, Phys. Lett. 62B, 253 (1976): Momentum distribution in the nucleus. 

[61] P. Jain, B. Pire and J. P. Ralston, Phys. Rept. 271, 67 (1996): Quantum color transparency and nuclear filtering. 

[62] A. S. Carroll et al. Phys. Rev. Lett. 61, 1698 (1988): Nuclear Transparency to Large-Angle pp Elastic Scattering. 



